Universality in open system entanglement dynamics 
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We show that the entanglement evolution of an open quantum system is the same for the vast 
majority of initial pure states, in the limit of large Hilbert space dimensions. 
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The ability to characterize and control quantum sys- 
tems with an ever increasing number of constituents is 
an ever more desired and necessary skill in many fields 
of modern physics [1], materials science [2], and even bi- 
ology [3]. From a strictly deterministic point of view, 
however, the rapidly growing generic complexity of such 
large systems apparently renders this task ineffable. Yet, 
similarly to thermodynamics, a statistical description of- 
ten allows for the extraction of robust, generic features 
which emerge in the limit of large system size, and imply 
quantitative predictions. 

Entanglement, an unmistakeable quantum signature, 
is a prime example of the above, apparent illusiveness of 
an exhaustive characterization as the system size is scaled 
up: The structure of many-particle entanglement turns 
more and more intricate with the exponentially increas- 
ing number of possible correlations between different sub- 
groups of particles. Thus, a complete characterization of 
a large, composite quantum state requires an experimen- 
tal overhead that increases exponentially with the num- 
ber of system constituents. Even worse, entanglement 
tends to get ever more fragile when enlarging the system 
size: the more degrees of freedom, the more difficult it be- 
comes to shield quantum coherences, which are necessary 
for entanglement, against the detrimental influence of a 
noisy environment. In such situations, the strong quan- 
tum correlations due to entanglement additionally need 
to be distinguished from classical correlations induced 
by the ambient noise. This is in general accomplished by 
high dimensional optimization procedures on the space 
of all quantum states [H [S] , leaving little hope for quan- 
titative predictions on entanglement evolution in large, 
and noisy systems. On the other hand, the signatures 
of entanglement in such systems are of high fundamen- 
tal and, potentially, practical interest, e.g. when it comes 
to harnessing the computational power of quantum algo- 
rithms [6 or assessing the role of quantum correlations 
in intrinsically noisy biological systems [3] . It is therefore 
a key issue and the subject of the present Letter to esti- 
mate the characteristic time scales in which entanglement 
is present in such adverse situations. 

Here we consider the fate of entangled states of large 
quantum systems, which require a high dimensional 
Hilbert space, in contact with an incoherent environment. 



We show that a statistical analysis over generic initial 
states unveils universal - state independent - open sys- 
tem entanglement evolution in the limit of large Hilbert 
space dimensions. In this thermodynamic limit, an effi- 
cient characterization of entanglement dynamics is thus 
again possible. 

More specifically, let us start with a composite quan- 
tum system in a pure state, characterized by a vector \x) 
in a Hilbert space H — Ha (Si Hb ■ ■ ■ (i> Hn with dimen- 
sion d = dAds ■ ■ ■ cLn- The system then undergoes some 
dynamics by virtue of its Hamiltonian, and also couples 
to uncontrollable degrees of freedom, which define its en- 
vironment (denoted by subscript E). When focusing on 
the system alone, its state p evolves as described by a 
time-dependent map A^: 

p(0) ^ p{t) = TvE [U{t)psEm{t)^] A,[p(0)], (1) 

where U{t) is the unitary time evolution operator ob- 
tained as solution of the Schrodinger equation for sys- 
tem and environment (subscript SE). Such maps Aj 
are only independent of the system's initial state p(0) 
if system and environment are initially uncorrelated or 
at most classically correlated [7j. This is the case for 
many implementations of quantum information tasks [8] , 
where initial system state are almost pure, and there- 
fore uncorrelated with the environment in very good ap- 
proximation. In what follows, we thus assume the ini- 
tial state of system and environment to be of the form 
Pse{^) = Ix)(xI'8)/3_e(0). Since, beyond this latter factor- 
ization, our ansatz makes no assumption on the specific 
form of the Hamiltonian and thus of U{t), and neither 
of /9£;(0), our subsequent results are applicable for arbi- 
trary open system dynamics of the given initial states, 
including non-Markovian effects that may arise in the 
course of the evolution. The following derivations are 
solely grounded in few geometrical properties of the set of 
states, and properties of entanglement measures thereon. 
The next steps and properties needed to prepare our main 
result are visualized in Fig. [l] 

As a measure for the relative effect of the open sys- 
tem dynamics Aj onto two states p and w, we choose the 
metric distance 

Dti{p,^) = ll/O-t^llxr := Tr|p-a;|, (2) 
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FIG. 1: State space is a convex set containing all quantum 
states of a given dimension. Open system dynamics due to a 
map At cannot increase the distance between any two states 
(as measured by the trace distance). Entanglement measures, 
as Ed, monotonically increase with the distance to the set of 
separable states (in gray). 

induced by the trace norm ||-||Tr on state space. The 
distance of two states measured by Dxr decreases mono- 
tonically upon application of any such Aj [5], 

i:»Tr[At(p),Af(w)] < 7?A, DTrlp,^) with ?7A, <1. (3) 

Note that, for closed dynamics, i.e., without coupling to 
the environment, equality holds with r^At = 1, as a con- 
sequence of the unitarity of quantum mechanics. Open 
dynamics, however, introduces additional mixing of the 
system state, such that the state space is effectively con- 
tracted. This dynamical effect is captured by the change 
of 77At . For example, when coupled to a thermal bath, all 
initial states converge to thermal equilibrium, and ?7Aj 
approaches zero asymptotically in time. 

Turning to entanglement, we do not focus on a specific 
entanglement quantifier, but rather only require that it 
slowly varies on the set of states. More specifically, we 
demand a "strong" form of continuity, namely Lipschitz 
continuity [10| . 

\E{p)- E{u)\<r^EDT,{p,Lo), (4) 

where jy^; is called the Lipschitz constant. Examples of 
such entanglement measures E[p) are the ones defined 
as the minimum distance of the state p to the set S 
of separable states: Ed{p) = mh\„(zs D{p,a) [5]. AH 
distances D which abide to the dynamic monotonicity 
condition ([3| [TT], fulfill the Lipschitz requirement with 
constant rjEo = li by virtue of the triangle inequality. 
Another example is negativity M , a computable entan- 
glement monotone |12| . but with a different constant "qj^f 
(see the appendix). 

Under conditions ([3| and Q, we can asses how much 
two initially pure states differ in their entanglement after 
exposure to the same incoherent dynamics. For states 
p{t) = At(|x)(xl) and uj{t) = At(|V)(Vl) we obtain 

\E{p) - E{uj)\ < ,7^i?Tr[A*(|x>(xl), At(|V')(VI)] 

<TiEm, D^,[\x){x\M){i^\] 

<27?i;77Aj|lx)-|V')||- (5) 



That is, the entanglement quantifier E inherits Lipschitz 
continuity, with constant 2riEr]At, even with respect to 
the Euclidean distance between the two initial state vec- 
tors in Ji. The difference in entanglement of the two 
mixed final states is essentially bounded by the distance 
between the initial states. With the knowledge of the 
entanglement of a single probe state, say of p(t), for an 
initial state |x), we are able to predict the entanglement 
evolution of any other pure initial state \ip) within an er- 
ror margin given by ([s]) - this without the need to evolve 
the state, and to compute the resulting mixed state's en- 
tanglement. 

Rather than studying specific instances of initial states, 
we aim at a statement about generic pure states in a sta- 
tistical manner. To avoid any bias, we use random states 
which uniformly cover the space of pure states. States 
that exhibit the properties of random states naturally 
appear in registers of quantum computers after long se- 
quences of gates |13j , in quantum systems with a chaotic 
classical counterpart [Hj, and also subsets of pure states, 
e.g. graph states which emerge in spin gases |15j . show 
the properties of a uniform distribution |16| . In addition, 
for many interesting scenarios of biology and chemistry 
such as transport phenomena in proteins, the initial state 
and/or the Hamiltonian are not exactly known, e.g. jl7j . 
In these cases a uniform distribution of initial states in 
the potential subspace is the best prior. 

Let us first note that the set of state vectors in H 
is isomorphic to a unit sphere S^''"^ in K^'*"^, where 
the real and imaginary part of the expansion coefficients 
into any basis constitute the coordinates, constrained by 
normalization and ignoring the global phase by choosing 
the first component real. With the above preparation, 
we can now infer our central result by employing the 
well-known Levi's lemma [THll20j : The probability for a 
deviation larger than e of E[p{t)] := E[At{\x) {x\)] from 
its mean {E){t) := / dtp ElAtllipjitpD] over aU initial 
states, given the uniform initial distribution on S^"*^^, 
exhibits an exponential suppression: (i) in the deviation 
e , and moreover, (ii) in the system dimension d: 



Pr {\E[p{t)] {E){t)\ >e)< 4exp f-C^^eA . 

(6) 

The constant C can be chosen (24 tt^)"^ [5T]. This holds 
for all entanglement quantifiers E that fulfill Q, even 
when quantifying multipartite entanglement. 

As a concrete application of this result, we consider 
the normalized negativity A/'/A/'max of a bipartite system 
on Ha ® Hb ■ Its Lipschitz constant (see appendix) reads 
^7 (A/'/A/'max) ~ dA/{dA — 1), and the following inequality 
for a deviation from the mean entanglement holds: 
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FIG. 2: (Color online) Distributions of negativity. The data for negativity is obtained for the least balanced partition (one vs. 

— 1 qubits) of = 3, 5, 8 qubits, at different snapshots of the evolution. We sample over 10 000 uniformly distributed pure 
initial states of A'^ qubits which evolve under local coupling to a decoherence reservoir. The system-environment interaction 
time t is parametrized by p = (1 — e^""*), with F the local decoherence rate. Clearly, the bigger A'^, the more concentrates the 
negativity distribution around its mean value, for all times, as exemplified in the logarithmically plotted standard deviation for 
system sizes A'^ = 2, . . . , 8 and three time steps of the evolution p = 0, 0.3, 0.5 (bottom right). Lines are linear fits to the data. 



Pr 



AA[p(t)]-(AA)(t) 



> e < 4cxp -C 



r 



,(2d^dB-l)(d^-l)2 



Ml 



(7) 



This expresses the concentration effect visuahzed in 
Fig. [2] for ensembles of — 2, ... ,8 qubits (two- level 
systems) , each locally coupled to an environment that de- 
stroys only the coherences. Negativity is here evaluated 
with respect to the least balanced bipartition of the qubit 
register, i.e., it quantifies the entanglement of one qubit 
{(Ia = 2) with the remaining TV — 1 qubits {ds — 2^~^), 
after application of the decoherence dynamics. We uni- 
formly sample 10 000 pure initial states, and parametrize 
the dynamics by the probability p for complete decoher- 
ence of a single qubit. Within the Markov approximation 
p ~ 1 — e^'"*, where F is the decoherence rate [55]. The 
log-plot of the variances for increasing system size, shown 
at the bottom right of Fig. [2] establishes the exponential 
concentration around the mean. Thereby, our numeri- 
cal results underpin the mixed state entanglement con- 
centration predicted by ([6]) and ([T]). Furthermore, the 
concentration observed by the numerical experiment is 
even stronger than estimated - the measured variances 
are approximately one order of magnitude smaller than 



the variance which we would infer from ([T]), with our 
above estimate of C and the Lipschitz constants. 

The above results determine the open system evolu- 
tion of entanglement with an error exponentially small 
in the dimension of the underlying Hilbert space (pro- 
vided r/E does not increase faster than Vd). Therefore, 
in high dimensions, it suffices to monitor the entangle- 
ment evolution of a single, generic pure state, in order 
to predict the fate of any other typical state subject to 
the same dynamics. Qualitatively different entanglement 
dynamics for different initial states ^23j will occur only 
as singular effects in sufficiently large quantum systems. 

Inasmuch as knowledge of a single state's entangle- 
ment evolution fully determines (in the present, asymp- 
totic sense) the result for arbitrary states, our result is 
reminiscent of previously derived entanglement evolution 
equations [24] l25], where the final entanglement of an 
arbitrary initial pure state was shown to be fully char- 
acterized by the entanglement evolution of a maximally 
entangled state. 
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The concentration of open system entanglement evo- 
lution, as spelled out by ^ and ([t]), has yet another 
bearing for the optimization problem routinely encoun- 
tered [HIS] when evaluating the entanglement of arbitrary 
mixed states. For mixed states generated through arbi- 
trary physical dynamics At from a uniform distribution 
of pure states, mixed state entanglement concentration 
suggests a reduction of the optimization space: A single 
representative of the obtained sample generated under a 
specific physical evolution^ selected with convenient prop- 
erties (symmetries in terms of its pure state decomposi- 
tion), will suffice to impose exponentially narrow con- 
straints on the optimization for all other states. 
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Appendix. In order to calculate the Lipschitz constant 
for the negativity of a bipartite system, we use its defini- 
tion through the trace norm of the partially transposed 
state [12, Af{p) := {\\(l(E) T){p)\\t, - l)/2, where T rep- 
resents the transposition acting just on Hb- By means 
of the triangle inequality, |||a;|| — \\y\\\ < \\x — y\\, and 
the linearity of the partial transpose we arrive at a first 
estimation for states p and w: 



value A/'max = (dA — l)/2, we finally obtain rn^js//js/^ 
dA/idA - 1). 
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Tr 



(8) 



The remaining norm can be estimated with the operator 
norm defined as ||A||op := sup^jdl Aa;||/||a;||). By choosing 
a particular x, not necessarily optimal, we obtain a lower 
bound: \\{l(g)T)\\op > r)(p-a;)||Tr/||p-a;||Tr ■ This 
leads to the next estimate for negativity: 

mp)~Miuj)\ < i||(i®r)||op||p-c.||T... (9) 

The maximization of the partial transposition is ob- 
tained with a maximally entangled state, which yields 
11(1 (8) T)||op = dA, with dA the dimension of the small- 
est subsystem of the bipartition. Consequently, we find 
Lipschitz continuity for negativity, 

\Af{p)^Af{0j)\<^\\p~0j\\Tr, (10) 
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